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Abstract. We derive new constraints on the neutron lifetime based on the recent Planck 2015 
observations of temperature and polarization anisotropies of the CMB. Under the assumption 
of standard Big Bang Nucleosynthesis, we show that Planck data constrains the neutron 
lifetime to r n = (907 ± 69) [s] at 68% c.l.. Moreover, by including the direct measurements 
of primordial Helium abundance of Aver et al. (2015) and Izotov et al. (2014), we show that 
cosmological data provide the stringent constraints r n = (875 ± 19) [s] and T n = (921 ±11) [s] 
respectively. The latter appears to be in tension with neutron lifetime value quoted by the 
Particle Data Group (r n = (880.3± 1.1) [s]). Future CMB surveys as COrE+, in combination 
with a weak lensing survey as EUCLID, could constrain the neutron lifetime up to a ~ 6 [s] 
precision. 
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1 Introduction 

The latest measurements of the Cosmic Microwave Background (CMB) anisotropies from the 
Planck satellite have presented a wonderful confirmation of the standard ACDM cosmologi¬ 
cal model [1, 2], The CMB is therefore providing a new “laboratory" where several aspects 
of fundamental physics can be tested and constrained. For example, the new Planck mea¬ 
surements have opened the interesting possibility of constraining the rate of the d(p, 7 ) 3 He 
reaction when combined with primordial deuterium measurements and assuming primordial 
standard Big Bang Nucleosynthesis (BBN) ([2, 3]). These cosmological constraints on nuclear 
reaction rates are comparable, in precision, with current direct experimental limits. Further¬ 
more, again thanks to the new precise data from Planck, the best constraints to date on 
atomic transition rates at last scattering have been obtained. Most notably, the amplitude 
of the 2s —> 1 two-photon rate that affects recombination is now determined by Planck data 
with a ~ 10% precision against the current ~ 43% experimental error (see [2] and references 
therein). 

In this paper we work along this line of investigation, showing that current cosmological 
measurements can provide interesting constraints on a key quantity in nuclear physics, i.e. 
the neutron lifetime. 

As it is well known, CMB physics is highly sensitive to the baryon density Ub and, to 
a lesser extent, to the primordial Helium fraction Y p . On the other hand, standard BBN 
can provide constraints on the primordial Helium fraction, once the baryon density and the 
neutron lifetime r n are fixed. At the same time, one can change perspective and use BBN to 
constrain the neutron lifetime, using the constraints on the primordial Helium fraction and 
the baryon density coming from CMB data. Clearly, this CMB+BBN constraint on r n can 
be further improved by considering the current estimates of Y p coming from astrophysical 
observations. 

The obtained constraints are, obviously, model-dependent. They relies not only on 
the assumption of standard BBN (see, e.g. [4]), but also on the several assumptions made 
under ACDM. Among them, probably the most important one is the assumption of a perfect 
knowledge of the amount of energy density in relativistic particles both at BBN and CMB 
epochs. This quantity, often parametrized by the number of relativistic degrees of freedom 
A^efj (see e.g. [2]), could be different from the standard value of N e g = 3.046 if sterile neutrinos, 
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axions or exotic light particles are present and could clearly affect the cosmological constraint 
on r n . 

On the other hand, this kind of analysis is rather timely since there is still no complete 
agreement on the experimental value of r n (for recent and comprehensive reviews see [5-7]). 
Indeed, while much technical progress has been made, measured lifetimes can vary by about 
a percent, depending on the experimental technique (see e.g. [5, 8]). Two main experimental 
approaches are available to measure the neutron lifetime: the “bottle method” and the “beam 
method”. In the “bottle method” a bottle is filled with Ultra Cold Neutrons (UCN), which 
are then emptied into a detector to count the number of remaining neutrons. Combining the 
five most recent UCN measurements [9-13] one obtains the very precise constraint of r^ >ottle = 
(879.6 ± 0.8) [s] (see [5]). On the other hand, in the “beam method” detectors record the f3- 
decay products in a part of the neutron beam. Combining the two most recent measurements 
[14, 15] based on the “beam method”, the following constraint is derived: r^ eam = (888.0 ± 
2.1) [s] (again, see [5]). As we can see, these two averages are in disagreement by ~ 8.5 s, 
i.e. by about ~ 3.7 standard deviations. Moreover, in over twenty years of experimental 
measurements of r n the average value has fluctuated over several standard deviations (for 
example, see discussion in [16]). The current experimental value quoted by the Particle Data 
Group (PDG hereafter) of r n = (880.3 ±1.1) [s] (see [16]) is an average over the best seven 
measurements, inflating the error by a scale factor of 1.9. 

It is therefore extremely interesting to investigate if current and future cosmological 
measurements could constraint the neutron lifetime to a precision that could be useful to 
shed light on present experimental discrepancies. At the same time, this analysis can show 
to what extent a precise experimental determination on T n could impact current and future 
cosmological constraints. 

Our work is organized as follows: in the next section we present our analysis method 
and we report the constraints on r n obtained from the Planck 2015 CMB data release. We 
also explore how current measurements of primordial Helium abundances can improve the 
bounds. In Section III we forecast the precision achievable by future CMB satellite missions 
while in Section IV we derive our conclusions. 

2 Current constraints on the neutron lifetime 
2.1 Method 

We consider the latest CMB anisotropy data coming from the Planck 2015 release (including 
polarization) and combinations of this dataset with CMB lensing ([17]) and Baryonic Acoustic 
Oscillation (BAO) surveys ([18]). Even if systematics may be clearly present (as indicated 
by the tension in the quoted values, see discussion below) we also consider the inclusion of 
recent astrophysical measurements of Y p . 

The analysis is based on the publicly available Monte Carlo Markov Chain package 
cosmomc [19] which relies on a convergence diagnostic based on the Gelman and Rubin statis¬ 
tic. We use the July 2015 version which includes the support for the Planck Likelihood Code 
2.0 [20] (see http://cosmologist.info/cosmomc/) and implements an efficient sampling of 
the space using the fast/slow parameters decorrelation [21]. We consider the ACDM model, 
we vary the six standard parameters: the baryon and cold dark matter densities Ub and w c , the 
ratio of the sound horizon to the angular diameter distance at decoupling 0, the reionization 
optical depth r, the scalar spectral index ng, and the overall normalization of the spectrum 
As at k = 0.05 Mpc -1 . We assume adiabatic initial conditions and we impose spatial flatness. 
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In our baseline analysis we assume standard radiation content (i.e. N e ff = 3.046), but we also 
consider the possibility of relaxing this assumption. As extra parameter, we consider the neu¬ 
tron lifetime T n . By itself, CMB is insensitive to r n but it depends on the primordial Helium 
abundance Y v (see e.g. [22]). For each step in the chain, we compute a primordial abundance 
y^BBN( w ^ Tn) obtained using a fitting formula based on the PArthENoPE BBN code [23]. Fix¬ 
ing the Helium abundance to Yp BBN therefore makes the CMB angular spectra sensitive to a 
variation in r n as we can see in Figure 1. 
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Figure 1 . Dependence of CMB power spectra on different values of r n , over plotted with the Planck 2015 
spectra. 

We present constraints on r n from cosmological data only and also in combination with 
4 He direct astrophysical observations. In the last decades many authors [25-32] have at¬ 
tempted to extrapolate primordial abundance of 4 He mainly from the observations of He I 
emission lines from low-metallicity HI I regions. For our analysis, we select results from Izotov 
et al. 2013 [29], Aver et al. 2013 [30], Izotov et al. 2014 [31] and Aver et al. 2015 [32], The 
first two papers represent the most recent results on Y p from the analysis of optical emission 
lines, updated with recent set of He I emissivities tabulated by [33] and larger sample of 
observed spectra from HII regions. Results from Izotov et al. 2014 and Aver et al. 2015 are 
the first results with a pioneeristic use of the NIR He I line at A = 10830A. For detailed 
description of the extrapolation of Y p and how they improve previous works see [29-32], 

Of course, a key point in our analysis is the precision to which primordial abundance of 
Y p is obtained assuming SBBN and given the values of Wf,, N e g and r n . In this work we rely 
on the PArthENoPE BBN code, that has been also used by Planck collaboration. However, we 
need to stress that other BBN codes are available and they may give slightly different values 
of primordial abundances. A recent comparison between available BBN codes has been made 
in [24], Comparing the value of Y p quoted in Tab. II of [24] with the version of PArthENoPE 
used in this analysis, we found a AYp = 0.0003, that corresponds to an accuracy of about 
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A r n ~ 1.5 [s]. As we will show later, this level of accuracy is acceptable when considering 
current constraints but it will clearly be an issue with future cosmological data. A second 
important point is the precision of the BBN fitting formula assumed in our analysis. In prin¬ 
ciple, for each values of Ub, N e g- and T n one should run a BBN code to compute the primordial 
abundances of light elements. However, a single run at fixed parameters is computationally 
rather expensive and it would be unrealistic to integrate it in the MCMC parameter explo¬ 
ration. For this reason, in BBN calculations, polynomial fitting formulas are used (for details 
see [2, 4]). Since the formula used in this paper is valid for values of r n close to the PDG one 
[16], we checked its precision for values within Ar n = 100[s]. Given a perfect knowledge of 
Ub and Ip, the expected precision on T n is 0.1 [s] for A T n = 20[s], 1.3[s] for Ar n = 60[s] and 
3.0[s] for A T n = 100[s]. We can then infer that these errors are negligible with respect to the 
statistical ones. 


2.2 Results from cosmological and astrophysical data assuming lV e ff = 3.046. 

We first report our constraints assuming a standard energy content in relativistic particles, 
i.e. N e ff = 3.046. The results of our analysis are reported in Table 1 where we report not only 
the constraints on T n but also, for completeness, the constraints on 1^ BN . We indeed point 
out that these constraints are obtained assuming a complete ignorance on the experimental 
value of the neutron lifetime. As we can see from Table 1, combining the low —£ temperature 
and polarization likelihood [20] with the small-scale temperature only likelihood, we obtain 
the constraint of r n = (918 ± 105) [s]. Including the small-scale polarization data from Planck 
provides a significant improvement, with r n = (907 ± 69) [s[. Further including to the CMB 
data BAO observations [18] we obtain r n = (915 ± 63) [s], while including also the lensing 
likelihood [17] we find r n = (894 ± 63) [s]. 


Dataset 

yBBN 

Ai[s] 

Planck TT 

0.254 ±0.021 

918±105 

Planck TT,TE,EE 

0.252 ±0.014 

907 ± 69 

Planck TT,TE,EE + BAO 

0.254 ±0.013 

915 ± 63 

Planck TT,TE,EE ± BAO + lensing 

0.249 ±0.013 

894 ± 63 


Table 1 . Constraints on YjP BN and T n at 68% c.l. based on the Planck 2015 release (see text). 


The constraints on r n coming from current cosmological data are in agreement with 
the PDG value but between one and two orders of magnitude weaker than those obtained 
by laboratory experiments. While future CMB measurements, as we discuss in the next 
section, could significantly improve these constraints, is also important to stress that the 
actual constraints are anyway obtained on very different scales in space and time. Assuming 
as current value for the neutron lifetime the PDG value of r BDG = (880.3 ± 1.1) [s] we see 
that the Planck data alone constrain any variation in T n from this quantity to be smaller than 
11% at 68% c.l. over a timespan of more than 13.5 Gyrs. 

We now provide cosmological constraints on the neutron lifetime by combining the lat¬ 
est CMB observations from Planck (including polarization) with four current astrophysical 
estimates of primordial Y p described previously. We combine those values with the Planck 
full dataset treating them as a gaussian prior on the input Helium abundance and running 
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the MCMC code as described above. We refer to the experimental direct measurement of 
primordial Helium as y^ ata . 

In Table 2 we show the results of combining the CMB data with the astrophysical 
determinations of Y p . We report the experimental value of Helium considered (Y^ ata ), the 
Helium mass fraction recovered from the MCMC chain (y^ BN ) and the values for the neutron 
lifetime (r n ) for all the datasets analyzed. As expected, the constraints on Y^ BN are now 
completely determined by the measured value of Y^ ata . 

The combined constraining power of the CMB data, sensitive to the baryon density, and the 
Helium astrophysical measurements allows to tight the bounds on r n by nearly one order of 
magnitude with respect to the CMB only constraint. 


Dataset 

■ydata 

P 

yBBN 

Ails] 

Planck TT,TE,EE + Izotov et al. (2013) 

0.254 ± 0.003 

0.2539 ±0.0029 

916 ± 15 

Planck TT.TE.EE + Aver et al. (2013) 

0.2465 ± 0.0097 

0.2484 ±0.0079 

888 ± 39 

Planck TT,TE,EE + Izotov et al. (2014) 

0.2551 ± 0.0022 

0.2550 ± 0.0022 

921 ± 11 

Planck TT.TE.EE + Aver et al. (2015) 

0.2449 ± 0.0040 

0.2455 ± 0.0038 

875 ± 19 


Table 2. Experimental priors lp data from different datasets of direct astrophysical observations. For each 
datasets we report the obtained value of Helium mass fraction yjf BN and the neutron decay time r„. 


Planck+Aver2015 

Planck+lzotov2014 

Planck+Aver2013 

Planck+lzotov2013 
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theoretical value 
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Figure 2. Values of r n , all reported with la error bars. In magenta we report the theoretical value. In 
green we report the combined values for the “bottle method” and “beam method” [5], the vertical violet stripe 
is the combined value [16]. In blue we report all values found including in the analysis the astrophysical 
measurements. 
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In Figure 2 we show all the values for the neutron lifetime evaluated considering individ¬ 
ually the astrophysical measurements. We also show in the same figure the combined values 
for the “bottle method” and “beam method” [5] , which are all used in the value quoted by the 
PDG [16], i.e. T n = (880.3 ± 1.1) [s], reported in the plot as well. 

Focusing on most recent astrophysical observations, we note that with Izotov et al. 
(2014) we obtain a value of r n in tension with the PDG one (about 3.7 a), while with the 
very recent Aver et al. (2015) measure we get r n = (875 ± 19) [s], in perfect agreement with 
the PDG but still not sensitive enough to discriminate between “bottle method” and “beam 
method”. The Izotov et al. (2014) and Aver et al. (2015) datasets provide values of r n , when 
combined with CMB data, that are in tension at the level of 2.1 a. This is expected since it 
reflects the tension between the two primordial abundance extrapolations. 

2.3 Results from cosmological and astrophysical data obtained varying A^ff- 

As mentioned above, we now check how our constraints are dependent on the chosen assump¬ 
tions for the dark radiation content. Therefore we include N e g- as extra parameter, performing 
the analysis both with cosmological data alone and in combination with primordial Helium 
abundance observations. Moreover it would be interesting to investigate if the inclusion of 
A r e ff could alleviate the tension on the value of r n between Izotov et al. (2014) and Aver et 
al. (2015). 

The results of this analysis are shown in Table 3 and in Figure 3. As expected we see an 
overall relaxation of the constraints on r„ (a factor of about 1.3 for the cosmological constraints 
and about 1.5 - 4 - 2 for the combined CMB and astrophysical measurements analysis), due to 
the well known degeneracy between A^g- and Y p . Given this degeneracy and considering that 
cosmological data now prefer values of A r e ff lower than the standard one, we obtain higher 
primordial Helium abundances and, as direct consequence, higher values of r n . Despite this 
overall behaviour, these new results are in agreement within 1 a with those discussed above 
(obtained fixing N e s = 3.046), as shown in Figure 3. In addition we always obtain values of 
A r e ff in agreement with the standard one within 1 a. As we can see, the tension between Aver 
et al. (2015) and Izotov et al. (2014) persists, albeit now at the level of 1.5 a. Of course, 
part of the previous 2 a discrepancy has been absorbed by the N e g parameter, resulting in a 
~ 0.5a shift between Izotov et al. (2014) and Aver et al. (2015) on N e g (see Table 3). 


Dataset 

yBBN 

Aeff 

7n[s] 

Planck TT,TE,EE 

0.263 ±0.018 

2.76 ±0.30 

986 ±109 

Planck TT,TE,EE + BAO + lensing 

0.258 ±0.014 

2.83 ±0.25 

953 ± 84 

Planck TT,TE,EE + Izotov et al. (2013) 

0.2542 ± 0.0029 

2.84 ±0.23 

933 ± 24 

Planck TT,TE,EE + Aver et al. (2013) 

0.2505 ± 0.0085 

2.90 ±0.25 

911 ±52 

Planck TT,TE,EE + Izotov et al. (2014) 

0.2552 ± 0.0022 

2.83 ±0.23 

939 ± 22 

Planck TT,TE,EE ± Aver et al. (2015) 

0.2458 ± 0.0040 

2.95 ±0.24 

884 ± 28 


Table 3. Constraints at 68% c.l. on neutron lifetime obtained varying N e a in addition to the other cosmo¬ 
logical parameters. In the last two rows we combine astrophyisical observations with Planck TT,TE,EE. 
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Planck+lzotov2014 


Planck+Aver2013 
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Planck+BAO+lensing 
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Figure 3. Values of r n , all reported with 1<t error bars. The vertical violet stripe is the combined value 
[16]. In green we report values obtained using cosmological data and in blue results from the combination 
of Planck TT,TE,EE with astrophysical observations of primordial Helium abundance. For comparison, we 
report results obtained fixing N e ft = 3.046 (dotted lines) and relaxing this assumption (solid lines). 

2.4 Theoretical estimate 

For reference, we also show in Figure 2 a theoretical estimate for r n = (883 ± 3) [s]. This 
value is the result including QED outer corrections (treating nucleons as a whole), and inner 
corrections depending on nucleon structure, the latter promoted at leading log level as in 
[34], Perturbative QCD corrections, weak magnetism and finite nucleon mass effects, related 
to proton recoil, as well as the Fermi function describing the electron rescattering in the 
proton Coulomb held, are also included. For more details on these points see e.g. [34-37]. 
The value of parameters entering this expression, such as vector and axial coupling Cy and 
Ca, fine structure constant, particle masses, strong coupling constant etc. has been updated 
to the most recent results reported in [16]. Finally, the low energy cutoff applied to the 
short-distance part of the 7 — W box diagram is chosen in the range Ma = 1 — 1.5 GeV. The 
theoretical error is obtained by propagating the uncertainties on all parameters to the value 
of r n , and is largely dominated by the present error on Ca/Cv. Of course, the theoretical 
estimate is still affected by smaller effects, such as for example, higher order terms in a, 
sub-leading corrections, residual average proton polarization due to parity non-conservation 
and so on [35]. All these further corrections are expected to be quite sub-dominant [38]. 
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3 Forecasts on future data 


Considering that current cosmological and astrophysical observations are not competitive with 
particle physics experiments yet and they also could be affected by unknown systematics, we 
now perform forecasts on forthcoming observations. 

First af all, it is interesting to evaluate the precision on Helium abundance from future 
astrophysical observations, combined with Planck, needed to reach the PDG accuracy. We 
find that to achieve 1 second precision on r n we need to measure Y v with an error of 0.0002, 
roughly an order of magnitude better than the current measurements. This would need 
a better control of the systematic errors on direct astrophysical measurements. Moreover, 
as discussed in the previous section, this would also need a significant, but in our opinion 
achievable, improvement in the numerical accuracy of current BBN codes, see e.g. a recent 
paper on possible corrections due to a better treatment of neutrino energy transport [39]. 


Experiment 

/sky 

Channel 

FWHM 

T /i A'• arcmin 

Q/U nK• arcmin 

COrE 

0.80 

105 

10’ 

2.68 

4.63 



135 

7.8’ 

2.63 

4.55 



165 

6.4’ 

2.67 

4.61 



195 

5.4’ 

2.63 

4.54 



225 

4.7’ 

2.64 

4.57 

AdvACT 

0.50 

90 

2.2’ 

7.8 

10.9 



150 

1.3’ 

6.9 

9.7 



230 

0.9’ 

25 

35 

SPT-3G 

0.06 

95 

1.6 ’ 

4.2 

5.9 



150 

1.0’ 

2.5 

3.5 



220 

0.68’ 

4.2 

5.9 

CMB-S4 

0.50 

150 

1.3’ 

1 

1.4 

CVL 

1.00 

150 

5’ 

0 

0 


Table 4. Experimental specifications for COrE [40], AdvACT [41], CMB-S4 [42], SPT-3G [43] and CVL. 
We report the used sky fraction / a k y , the used channels, the beam FWHM in arc-minutes and the sensitivity 
per pixel in temperature and polarization. 


It is also interesting to evaluate the future constraints on r n achievable from an anal¬ 
ysis based just on cosmological data as CMB or galaxy clustering. CMB sensitivity on r n 
is encoded in the small-scale region of the power spectrum, therefore we can expect tighter 
constraints from next generation CMB missions, planned to measure such ^ range. To fore¬ 
cast the sensitivity on r n that future CMB experiments can achieve, we simulate different 
synthetic datasets. We consider the following future experiments: COrE [40], AdvACT [41], 
Stage IV CMB experiment (CMB-S4 hereafter) [42] and SPT-3G [43]. All the experimental 
specifications are reported in Tab. 4. Concerning the COrE experiment, we have chosen the 
5 channels in the range 100-220 GHz. We also consider an ideal Cosmic Variance Limited 
(CVL hereafter) case, to evaluate the most accurate precision that can be reached with CMB 
experiments, assuming no systematic effects and variance dominated by the signal in the 
considered range of multipoles. In order to extract the cosmological parameter constraints 
we generate simulated datasets following the, now common, approach described in [44]. We 
assume, as fiducial model, the Planck 2015 best fit. Moreover we assume that uncertainties 
due to foreground removal are smaller than statistical errors, that beam uncertainties are 
negligible and white noise. As Likelihood function we use the Inverse Wishart distribution, 
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reducing the degrees of freedom at each £ by f s k y and ignoring the correlations. Results are 
reported in Table 5. 


Dataset 

yBBN 

T n[s] 

Planck TT,TE,EE + AdvACT 

0.2464 ± 0.0065 

879 ± 32 

Planck TT.TE,EE + CMB-S4 

0.2475 ± 0.0037 

884 ± 18 

Planck TT,TE,EE + SPT-3G 

0.2487 ±0.0091 

890 ± 44 

COrE 

0.2467 ±0.0023 

880 ± 11 

CVL 

0.2467 ±0.0011 

880.7 ± 5.5 

Planck TT.TE,EE + Euclid 

0.2521 ± 0.0069 

907 ± 34 

COrE + Euclid 

0.2467 ±0.0014 

880.3 ±6.7 


Table 5. Future constraints on y p BBN and r n at 68% c.l. achievable from the considered cosmological datasets 
(see text). 


Combining the Planck 2015 likelihood with simulated datasets for the future ground 
based missions SPT-3G, AdvACT and CMB-S4 the 1 a error on the measured r n is reduced 
respectively to 44 [s], 32 [s] and to 18 [s]. The next generation satellite mission COrE could 
lower this constrain to 11 [s] while the maximum intrinsic limit achievable with the CMB data 
alone is about 5.5 [s]. Finally, we also perform a forecast for the future Euclid mission [45], 
combining it with the latest Planck data and with the future COrE mission. We use the 
Fisher matrix formalism, following [46] and the experimental specifications reported in [47]. 
We sum the Euclid Fisher matrix with the CMB inverse covariance produced by cosmomc, 
assuming gaussian probability distribution for the cosmological parameters. Despite Euclid 
being not sensitive to the variation of r n , it clearly helps in breaking the degeneracies between 
the cosmological parameters. Infact we obtain a significative reduction of the error on Yp BBN 
and T n with respect to the Planck only and CorE only cases, as reported in Table 1. 

4 Conclusions 

In this paper we have presented new constraints on the neutron lifetime combining a set 
of current cosmological data. In particular, we have shown that the recent Planck 2015 
temperature and polarization data release, under the assumption of standard BBN, constrains 
the neutron lifetime to r n = (907 ±69) [s]. While the uncertainty on this value is significantly 
larger than current experimental estimates, it provides a new and completely independent 
measurement. 

We have combined the Planck constraints with current astrophysical determinations 
of Yp. We have found that in this case a A r n ~ 10 4- 20 [s] accuracy could be reached. 
This level of sensitivity is still not able to discriminate between “beam method” and “bottle 
method”. Anyway, forecasting future constraints achievable by a combination of Euclid and 
COrE± satellite mission, we have shown that these experiments could reach,in principle, a 
A t„ ~ 6 [s] sensitivity, providing also an improved accuracy in current BBN codes. Therefore, 
by combining future and more precise cosmological and astrophysical observations we could 
a reach sensitivity comparable with current experimental uncertainties, with the caveat that 
all the possible sources of systematics are properly taken into account. 
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